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AN ACTUATOR EXTENSION TRANSFORMATION FOR A 
MOTION SIMULATOR AND AN INVERSE TRANSFORMATION 
APPLYING NEWTON-RAPHSON'S METHOD 

By James E. Dieudonne, Russell V. Parrish, 
and Richard E. Bardusch 
Langley Research Center 

SUMMARY 

A set of equations which transform position and angular orientation of the centroid 
of the payload platform of the six -degree -of ^ -freedom motion simulator at the Langley 
Research Center into extensions of the simulator's actuators has been derived and is 
based on a geometrical representation of the system. An iterative scheme, Newton - 
Raphson’s method, has been successfully used in a real-time environment in the calcu- 
lation of the position and angular orientation of the centroid of the payload platform when 
the magnitude of the actuator extensions is known. Sufficient accuracy is obtained by 
using only one Newton -Raphson iteration per integration step of the real-time environment. 

INTRODUCTION 

To enhance the capability of producing realistic aircraft simulations, in December 
1971, NASA Langley Research Center acquired a Singer Simulation Products Division 
six -degree -of -freedom motion simulator. This paper describes two problems inherent 
to the design of this particular motion simulator and the methods used to solve these 
problems. First, since the motion simulator does not have independent drive systems 
for each degree of freedom but achieves motion in all degrees of freedom by a combina- 
tion of actuator extensions, the base cannot be driven with the usual inputs of position 
and angular orientation of the centroid of the payload platform. This fact requires that 
a transformation be developed which will convert motion cues into actuator extensions. 
This transformation will be discussed first in this paper. 

For the purpose of evaluating the performance of the simulator hardware, an inverse 
' transformation is desirable. The actual position and angular orientation of the centroid of 
the payload platform would be obtained from the transformation of the actual magnitudes 
of the actuators. A comparison of desired position and orientation of the centroid with 
actual position and orientation would then yield the error in the base drive and provide a 
method of determining dynamic servo performance. This inverse transformation is also 



required for optimal washout filter design; therefore, the development of a method of 
computing this inverse transformation that will operate in a real-time environment is 
presented. 


SYMBOLS 


Values are given in both SI and U.S. Customary Units. The measurements and 
calculations were made in U.S. Customary Units. 


A i 


A; 


i,m 


vector in fixed coordinate system from origin of moving coordinate system 
to upper attachment point of actuator i 

known vector in moving coordinate system from origin of moving coor- 
dinate system to upper attachment point of actuator i 


a i,mi’ a i,ni2’ a i,ni3 


elements of Ai 


Bi 


known vector in fixed coordinate system from origin of fixed coordinate 
system to lower attachment point of actuator i 

elements of B; 


b; ,b; ,bi 

H x 2 *3 
Ej extension of actuator i 


f( ) 


H 


vector function 


height above lower bearing plane 


i m ,jm,k m components of unit vectors defined in moving coordinate system 


i 0 ,j 0 ,k 0 components of unit vectors defined in fixed coordinate system 
*i 


vector in fixed coordinate system from lower attachment point to upper 
attachment point of actuator i 


c. c ■ C- 
i,x» i,y» i,z 


components of 


*ii 


actual magnitude of actuator i obtained from instrumentation 
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Ip- I magnitude of actuator i when payload platform is in its neutral position 

I 1 1 neut 

O null vector 

R vector in fixed coordinate system from origin of fixed coordinate system 

to origin of moving coordinate system 

jk vector in fixed coordinate system from origin of fixed coordinate system 

to upper attachment point of actuator i 

[t] Euler angle of transformation matrix 

Tjj elements of T matrix 

x,y,z elements of R which determine the inertial position of payload platform 

~a unknown parameter vector which is vector root of f( ) 

\p,6,(p inertial angular orientation of payload platform 

Subscripts: 

m moving coordinate system 

max maximum 

min minimum 

n past value 

n+1 new value 

neut neutral 

Superscript: 

T transpose 

A bar over a symbol denotes a vector. 
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MOTION BASE DESCRIPTION 

The motion base (fig. 1) consists of a payload platform upon which a cockpit and 
visual display will be mounted. The payload platform is driven by six hydraulically 
powered position servomechanisms forming a synergistic six-degree -of -freedom motion 
system. As shown in figure 2, the points where the actuators connect to the payload plat- 
form form an equilateral triangle approximately 3.66 m (144 in.) on a side. The actua- 
tors have a minimum length of 2.62 m (103 in.) and a maximum length of 4.14 m (163 in.) 
giving a fully settled height of the platform of approximately 2.05 m (81 in.) and a fully 
raised height of 3.99 m (157 in.) measured from the floor. The system design allows 
for a 6804 kg (15 000 lb) payload and provides the limits of performance shown in table I 
from the neutral position. 

ACTUATOR EXTENSION TRANSFORMATION 

Because of the design concept of the hardware, activation of all six hydraulic actu- 
ators is generally required for motion in each degree of freedom. Thus, the simulator 
must be driven with actuator extensions, but since aircraft simulations provide position 
and angular orientation of the centroid of the payload platform, a transformation to actu- 
ator extensions is necessary. 

Two coordinate systems are defined as shown in figure 3; one has its origin at the 
centroid of the lower fixed platform of the simulator and the other, at the centroid of the 
moving payload platform. The coordinate systems are both right-hand systems and their 
axes are alined when the payload platform is at its neutral position. When one actuator 
at a time is considered, figure 4 defines the vector relationships between the origins of 
each coordinate system and the actuator attachment point of each platform. Figure 5 
shows these relationships for a general orientation of the motion simulator. These rela- 
tionships yield the vector equations 

(1) 

which after subtraction become 

= Aj + R - Bj (2) 

This equation is defined with respect to the fixed reference frame. Since the known vec- 
tors A^ m are defined with respect to the moving frame, one must apply an Euler angle 
transformation in order to determine the corresponding vectors Aj in the fixed refer- 
ence frame. By using a i//,0,0 order of rotation, the transformation is 


r j = A^ + R 
h = Bi + li 


4 


where 


im 


V 

3m 

s 

II 

3o 

^m 


k o 



T 11 T 12 

T 13 
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s 

t 21 t 22 

T 23 

i 




[_ t 31 t 32 

t 33 



and 





Tj j = cos i p cos 8 

T 12 

= sin \p cos 0 

Tjg = -sin 0 

T21 - cos \p sin 8 sin <p 

t 22 

= sin \p sin 0 sin 0 

T23 = cos 0 sin 0 

- sin ip cos 0 


+ COS 1 p COS 0 


T31 = cos \p sin 8 cos 0 

t 32 

= sin 10 sin 0 cos 0 

T33 = cos 0 COS 0 


+ sin \p sin 0 - cos 0 sin 0 


Applying this transformation to the Aj^ m vectors yields 

Si=M T A i>m 

Substituting this relation into equation ( 2 ) yields 

^i = [ T ] TX i,m+R -Bi (3) 

which can now be solved since m and Bj are known constant vectors, and R and 

M T are calculated from the given values of x, y, z, ip, 8 , and 0 . The length of 
each actuator is then 
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and the actuator extensions Ej are defined by 



where I h I is the known neutral position value of actuator i. 
I 1 meut 

INVERSE TRANSFORMATION 


The actuator extension transformation permits the driving of the motion simulator 
with position and angular orientation inputs, but a means for determining the simulator's 
response to these signals and its positioning accuracy was unavailable. A solution would 
be to take the actual actuator extensions, available from potentiometers mounted on each 
actuator, and transform them into the position and angular orientation of the centroid of 
the payload platform. Actual error, phase lag, and so forth, could then be determined 
from comparisons of the commanded and actual position and orientation. This transfor- 
mation is simply the inverse of the actuator extension transformation. This inverse 
could be obtained if the vectors were available, but the potentiometer reading yields 
only the magnitude of the corresponding actuator [ j and not the required vector. 
With only this information available, the problem becomes that of solving six simulta- 
neous nonlinear equations for six unknowns (x, y, z, if/, 6, and <p ). 


The approach was to apply an iterative numerical method known as the Newton- 
Raphson method. (See ref. 1 (pp. 447-453) and ref. 2 (pp. 19-26).) For a vector-matrix 
equation, this method is a general method of computing the vector root of 
equation (4) 

f(« ) = O (4) 

The iteration formula has the form 


Q! 


n+1 


= a n - 


^Qn) 

da n 


-1 

f "K) 


(5) 


In order to apply this method to the problem, a function satisfying equation (4) must first 
be defined. Since the actual magnitudes of the actuators are known, a function fj(a) can 
be defined as 



( 6 ) 
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where 


x 


y 


a = 


z 


e 


<t> 


f i (°0 

f 2(“) 



and [ -Cj t a is the actual magnitude of actuator i obtained from instrumentation. Sub- 
stitution of equation (3) into equation (6) and performing the vector multiplication yields 
equation (7). Note that each vector multiplication results in a scalar and permits the 
grouping of terms as follows: 



T rp 2 

+ R R + B t Bi - | i i | a (7) 
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By expanding in terms of elements, equation (7) becomes 


f i(“ ) - »!%,, + a f,m 2 + a i,m 3 + b , 2 l + b ? 2 + bf 3 + x 2 + y 2 + z 2 - | tj \\ 


* 2 ( X ' bi l)( ai . m l Tl1 + a i,"> 2 T 21 + a i,m 3 T 3l) 
+ 2 ( y ' b ‘2)( ai > m l Tl2 + ill . nl 2 T22 + a i,nl 3 T 32) 
+ 2 ( z ' b i 3 )( a i.m 1 T 13 + a i,m 2 T 23 + a t,m 3 T 33) 

- K x \ + y \ + z \) 


( 8 ) 


Taking the partial derivatives of f^a) with respect to the elements of a yields 


«i(«) 


3x 


af i( a ) 




9f i ( Q! ) 


dZ 


fiW 


= 2 ( x + a i,m 1 T H + a i,m 2 T 21 + a i,m 3 T 31 ~ 

= 2 (j + a i,m 1 T 12 + a i,m 2 T 22 + a i,m 3 T 32 ' b i 2 ) 

= 2 ( z + a i,m 1 T 13 + a i,m 2 T 23 + a i,m 3 T 33 “ h i 3 ) 

= ' 2 ( X " + a i,m 2 T 22 + a i,m 3 T 32) 

+ ' bi 2)( ai > m l Tl1 + a i> m 2 T21 + ai ’ m 3 T31 ) 


(9a) 


(9b) 


(9c) 


(9d) 
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The solution to the problem is now complete, since for a given set of initial conditions of 
a and corresponding solutions to equations (8) and (9), equation (5) can be solved. The 
final form of the solution is 
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where Q J is defined as a function of past value of oi, Qf n ; £ Jn+1 * S ^ e ^ nec * as 
new value of a, « n+1 ; and f t implies fi(a). With an appropriate set of initial con- 
ditions, equation (10) is repeated until some predetermined convergence criterion (desired 
accuracy) is met. 

The fact that the inverse transformation must operate efficiently in a real-time 
environment places additional constraints on the solution of the inverse problem. Besides 
assisting in the validation of hardware performance, the inverse transformation is to be 
used in the feedback loop for optimal washout techniques. Figure 6 illustrates the major 
parts of a real-time motion simulation by using a washout technique requiring feedback. 

As shown in the figure, the motion simulation requires the addition of all the starred 
blocks to a normal fixed-base simulation. These additions will increase the computational 
time required per integration step for all motion simulations. The inverse transforma- 
tion, being an iterative method, requires a variable amount of computational time based 
on the number of Newton-Raphson iterations; therefore, the minimum number of iterations 
required to give the desired accuracy had to be determined. 

RESULTS 

Actuator Extension Transformation 

The actuator extension transformation was shown to be valid by the use of several 
independent checks. 


Inverse Transformation 

To operate efficiently in a real-time digital environment, the inverse transformation 
must yield the desired accuracy in the minimum amount of computational time possible. 
Since most of the real-time digital simulation programs at the Langley Research Center 
utilize an integration step size of 1/32 second, the inputs to the inverse transformation 
| | will be changing every 1/32 second. Based on the maximum servo drive rates of 

the motion base and the sampling rate of the actuator magnitudes, the minimum number 
of Newton-Raphson iterations per sample required to give the desired accuracy had to be 
determined. The desired accuracy was defined to be anything within the position specifi- 
cations of the servo drives of the motion base. 

The data necessary to determine the minimum number of Newton-Raphson itera- 
^ Iqus acceptable were obtained in the following manner* The actuator extension transfor- 
mation was initialized at a particular starting point and an input to one degree of freedom 
was swept through the allowable range of that degree of freedom with the maximum allow- 
able rate (determined from the maximum servo drive rates). The resulting actuator 
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extensions were used to drive the inverse transformation; thus, a motion base with perfect 
response was simulated. Any differences that existed between the inputs to the actuator 
extension transformation and the outputs of the inverse transformation would then be 
attributable to the number of Newton -Raphson iterations. 

Although many different starting points were tried, no attempt was made to conduct 
a Monte Carlo simulation of the infinite number of starting points available. The results 
to be presented represent the worst case encountered. 

Table II presents the maximum error obtained in any degree of freedom by using 
one, two, and three iterations for an input in each individual degree of freedom. For 
example, with only the x-input changing (at its maximum allowable rate), the maximum 
error that occurred during the sweep of the x-range was in z, 6.667 x 10"® m 
(2.624 x 10“ 3 in.) for one iteration. The results demonstrate that with one iteration 
the maximum errors obtained, 9.287 X 10"** m (3.655 X 10“^ in.) in translation and 
2.719 x 10-5 rad (1.558 x 10 3 deg) in rotation, are within the specifications of the servo 
drives, 3.074 x 10" 3 m (0.121 in.) in translation and 8.726 x 10“^ rad (0.05°) in rotation. 

Table III shows the maximum error obtained in each degree of freedom, by using one 
iteration, for maximum rate inputs in individual degrees of freedom. Maximum errors 
for one-half and one-fourth maximum rates are also shown along with maximum error 
obtained from maximum rate in all six degrees of freedom simultaneously. 

It is interesting to note that the maximum error in the translational degrees of 
freedom, regardless of the input, always occurs in z. From a comparison of these runs, 
it can be determined that the error is directly proportional to rate of drive and that even 
the error incurred from the physically unobtainable case of maximum rate in all degrees 
of freedom simultaneously is still within the required servo drive specifications. 

It should be noted that the addition of the inverse transformation with one Newton- 
Raphson iteration to a real-time digital simulation program will increase the total com- 
putational time per integration step by 0.0018 second and the total memory by 1264 octal 
core locations. 


CONCLUDING REMARKS 

A transformation which provides actuator extensions for a six -degree -of -freedom 
motion base at the Langley Research Center from the motion cues of a real-time digital 
aircraft simulation has been developed and implemented. 

An inverse transformation which provides position and angular orientation of the 
centroid of the payload platform from the measured magnitude of the actuator extensions 
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has been developed by using a Newton -Raphson technique. The transformation has been 
shown to yield the required accuracy with one Newton-Raphson iteration for successful 
operation in a real-time environment, with only a slight increase in memory and com- 
putational time. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 19, 1972. 
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TABLE I.- PERFORMANCE LIMITS 


Degree 

of 

freedom 

Performance limits 

Position 

Velocity 

Acceleration 

Longitudinal (x) 

Forward 1.245 m 
Aft 1.219 m 

±0.610 m/sec 

±0.6g 

Lateral (y) 

Left 1.219 m 

Right 1.219 m 

±0.610 m/sec 

±0.6g 

Vertical (z) 

Up 0.991 m 

Down .762 m 

±0.610 m/sec 

±0.8g 

Yaw (0) 

±0.559 rad 

±0.262 rad/sec 

±0.873 rad/sec^ 

Pitch ( 0 ) 

±0.524 rad 
-.349 rad 

±0.262 rad/sec 

±0.873 rad/sec^ 

Roll (0) 

±0.384 rad 

±0.262 rad/sec 

±0.873 rad/sec^ 


TABLE II.- MAXIMUM ERROR OBTAINED FROM NEWTON -RAPHSON'S ALGORITHM 


Degree of 
freedom 
of input 

Maximum error after - 

1st iteration 
(*) 

2d iteration 
(*) 

3d iteration 

X 

(z) 6.667 X lO" 5 m 

(z) 8.160 X lO" 10 m 

(z) 2.925 X 10- 14 m 

y 

(z) 6.667 X 10” 5 m 

(z) 8.160 XIO-iOm 

(z) 2.415 X 10-1 4 m 

Z 

(z) 9.287 X 10' 5 m 

(z) 3.183 X lO' 10 m 

(z) 4.621 X 10 -14 m 

0 

(0) 1.926 x 10' 5 rad 

(0) 1.109 X lO' 10 rad 

( 4 /) 9.920 X 10“ 15 rad 

e 

(6>) 2.119 X 10‘ 5 rad 

( 6 ) 8.980 x 10" 11 rad 

(6) 1.984 x 10' 15 rad 

0 

(0) 2.719 X 10' 5 rad 

. 

(0) 9.245 x 10 -10 rad 

(0) 1.745 X 10- 17 rad 


* Quantity in parentheses denotes the degree of freedom in which maximum 


error occurred. 
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9.044 X 10" 6 7.550 X 10” 6 2.619 X 10" 5 8.841 X 10” b 1.178 x 10” b 1.205 x 10” 


























Figure 1.- Six-degree -of -freedom motion simulator. 
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Figure 4.- Vector relationships for actuator i 
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Figure 5.- Particular orientation of motion simulator. 



19 


Control Data 6600 



20 


NASA-Langley, 1972 8 


L-8505 


motion simulation. Starred blocks indicate additional software required. 





